On the kinetic theory of turbulence in plasmas 
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From the nonlinear Vlasov equation, a turbulence scattering term is found to 
describe the stachasticity on the time scale longer than the turbulence correlation 
time. The evolution of the plasma distribution is determined by the well-understood 
unperturbed motion of charged particles, with the effects of fluctuating part of fields 
described by the turbulence scattering term, from which, one can identify various 
important physics, covering from the linear and quasilinear to the nonlinear regime, 
in particular the collections between the widely used Kadomtsev-Pogutse's equation 
and Frieman-Chen's equation. 
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Plasma transport is one of the most challenging scientific problems. It is theoretically 
treated by solving the kinetic equation that determines the evolution of the particle distri- 
bution in the phase space, and the kinetic equation is nonlinearly coupled with Maxwell's 
equations, because the motion of charged particles depends on the electric magnetic fields 
and meanwhile the fields depend on the distribution of the charged particles [iKD]. Gener- 
ally, the fields can be decomposed into the fluctuating part and the averaged part jlj, and 
the motion of charged particles in averaged fields is usually well-understood. The fluctua- 
tion of fields can be due to the particle discreteness and the collective instabilities of the 



plasma, which are known as collisional process 
Macroscopically, both the collisional transport 
be described by diffusion-convection fluxes. 



and micro-turbulence js-2|, respectively. 
-|4j and the turbulent transport p-|9| can 



Up to the present, the collisional transport and the turbulent transport are theoretically 



investigated in different ways. In the well-established collisional transport theory 



4|, the 



particle motion due to the fields fluctuation induced by the particle discreteness is decoupled 
from the unperturbed motion in the averaged fields [l[ , and the decoupled motion is described 
by the Landau collision operator jsj], which can be written in a Fokker-Planck form. In the 
turbulent transport theory and numerical simulation, however, the particle motion due to 
the field fluctuation is usually not decoupled, and the particle trajectory has to be evaluated 
with the full fields js, 7|. In this way, one has to solve the nonlinear Vlasov equation or its 
gyro-averaged version for the magnetic confinement fusion plasma, by following the particle 
trajectory in the full fields (lo|4l3|] pj-|l8|. 



Recently, it has been proposed that the particle motion in the fluctuation fields due to 
the collective process can be decoupled from the unperturbed motion, by using the Hamil- 
tonian Lie-transform perturbation method, and the quasilinear transport theory has been 
successfully demonstrated {J. In this paper, by using this new method, we will show that, 
in a nonlinear turbulent plasma, the particle motion due to the fluctuating part of fields can 
be decoupled from the well-understood unperturbed motion, and the effects of fluctuating 
fields on the particle distribution can be put into a turbulence scattering term, which is 
written in a Fokker-Planck form as a divergence of diffusive- convective flux in the phase 
space. The turbulent transport and the collisional transport can be theoretically treated in 
a unified way. 



3 



We begin with the general kinetic theory of plasma transport 



a. 



V (/„) = (), (la) 



e 



V = d t + v ■ d x + — (E + v x B) ■ d v , (lb) 



V ■ E = 1 (^ Pex + S s J d?ve s f)j , (2a) 

V x B = no (j ex + S s y d 3 ve s vf s J + -^£7, (2b) 
Vx£= -a t f?, (2c) 

V • £? = 0, (2d) 

where f s (x, v, t) is the distribution function of particle species s at time t, with x and u the 
particle position and velocity, respectively; e s and m s the charge and mass of the particle, 
respectively; E and B the electric and magnetic fields, respectively. p ex and j ex are electrical 
charge density and current density due to the external sources, respectively. In the following, 
the subscript s shall be omitted for simplicity. 

For comparison, we write down the kinetic equation for collisional transport 



V (/) = C (/) = -d v 



ac - 7^d c ■ d v 



(3a) 



V = d t + v ■ d x + - (E + v x £? ) • d v , (3b) 
m 



where the velocity space convection vector ac and diffusion tensor dc depend on [3| /. Eq 
and B Q are averaged fields, with the fluctuation due to particle discreteness removed to 
the collisional scattering term, C (f). Usually, the effects of the micro-instabilities are not 
included in collisional transport theory. 

In the following, to concentrate on the turbulent transport, we shall deal with Eq. (fTaj) . 
with E = E + 5E and B = B + SB understood as the fields without fluctuation due to 
the particle discreteness, and Eq and Bo understood as the equilibrium fields that do not 
include the fluctuations (5E, SB) due to micro-instabilities. 

We begin with the particle motion in the unperturbed electromagnetic fields. The un- 



perturbed fundamental one-form (phase space Lagrangian) 19j is written in terms of the 
noncanonical variables z = (x, v) as 

7o = [mv + eA (x)] ■ dx — H (x, v) dt, (4) 



i 



with the unperturbed Hamiltonian H = (l/2)m\v\ 2 + e<fr (x). E = — V$ - -Bo = V x A . 
The well-unerstood equations of unperturbed motion are 

^ = {z i ,H } = ji i d j H , (5) 



with Jq = {z l , z^}, and the unperturbed Poisson bracket 



1 e 
{/, 9} = — (d x f ■ d v g - d v f ■ d x g) -B ■ d v f x d v g. (6) 



m ' w? 



Introduce the perturbations of the electric field SE = —V5§(x,t) — d t 5A(x,t) and the 
magnetic field 5B = V x 5 A, which are due to the micro-instabilities. The fundamental 
one- form is separated into the unperturbed part and the perturbation part, 7 = 70 + 71, 
with 

71 = e5A(x, t) ■ dx- e5&(x, t)dt. (7) 



We shall endeavor to decouple the perturbed motion due to the turbu 



unperturbed motion by using the recently proposed method of I-transform 9| , which, as a 



ence from the 



special Lie-transform 



20 



21] , makes the equations of motion in the perturbed fields formally 
identical to the unperturbed motion. The ordering parameter is e$ ~ | \ ~ | ^| <C 1, 
with Vth the thermal speed of the particle. To 0(e$), the phase-space Lie transform 20|, \21\ 
from the old variables z to the new variables z is 

£ = z l + G\{z) + G\ + -GidjGl (8a) 



^_ G i^_Gi + l G { djG l ( 8b ) 



where G\ and G2 are the first order and the second order Lie-transform generating vectors, 
respectively. 

The distribution function (a scalar) is transformed according to / (z) = f (z), which can 
be written up to 0(ef) as 

/ = / - (G\ + G\) dj + -GidiGifyf, (9a) 
/ = f + (G\ + Gi) dj + l -G%G{dJ. (9b) 

The straightforward calculation of the I-transform begins with 70 and 71, and it is similar 
to Ref. y. The results needed are summarized as follows. 
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The gauge functions of the I-transform, S n , are solved by 

V Q (S n ) = e(5V n -v5A n ). (10) 
For the low-frequency turbulence, the gyro-phase dependent part of S n can be analytically 



solved js] to decouple the fast gyro-motion. 

The I-transform generating vector can be written as 



Gl = --d v S n , (11a) 
m 

G v n = -d x S n + -^B x d v S n + -5An. (lib) 



m m 2 m 



The effective potentials are given by 

[5Ai,5V 1 ] = [5A,5$], (12a) 
[6Aa,6* 2 ] = [^Gf x 5B, l -G x x ■ SE}. (12b) 

With the I-transform generating vectors given above, the transformed fundamental one- 
form is identical to the unperturbed one, 

To = [mv + eA (x)] ■ dx — H (x, v) dt, (13) 

and hence the equations of motion are simply 

dz { /dt = {z\ H } = JtfdjHo. (14) 

The Vlasov equation is transformed to 

V (f) = 0. (15) 



It is not hard to show that G n 's are incompressible flows in the phase space 



a, 



d z ■ G n = 0. (16) 

To proceed, we make some remarks on the I-transform. G n 's describe the excursion 
from the unperturbed orbit. In describing the turbulence wave scattering, Eq. (1101) may 



be understood as the stochastic equations 



6j. It is important to understand that the I- 



transform is a perturbation theory; its validity depends on the condition that the real orbit 
does not deviate much from the unperturbed orbit. When the perturbation is weak, this 
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condition is well satisfied; clearly, this is the case of linear and quasi-linear theory discussed 
in Ref. 0. However, in the nonlinear stage, the perturbation of fields can be strong, and 
hence the real orbit of the particle may deviate significantly away from the unperturbed 
orbit in a long-time. Therefore, in the nonlinear stage, one can only apply the I-transform 
in a small time interval. 

To evolve the particle distribution within a small time interval At, we denote the dis- 
tribution function at time t — At as / (z,t — At). To apply the I-transform to find the 
distribution function at time t, one should keep in mind that it is the functions rather than 
the values that the Lie-transform operates on 20|, |21 |. 

Set S n (z, t — At) = 0, we have G n (z, t — At) = 0. This means that we do not make any 
transform at time t — At; therefore, according to Eq. ( |9a|) . we have 



f(z,t-At) = f(z,t-At) 



(17) 



Integrating Eq. (TTOj) along the unperturbed orbit from t — At to t, we can find S n (z, t) 
and G n (z,t). Since this is an infinitesimal I-transform, we denote the infinitesimal gauge 
functions and the infinitesimal generating vectors as s n (z,t) and g n (z,t), respectively. In- 
tegrating Eq. (EE5J), we found 



f(z,t) = f[z-Az(z),t-At] 

ft- At 

-Az(z) = , dt{z,H }. 



Subsitituting Eq. ( 1TT1) . we found 



f(z,t) = f[z-Az(z),t-At] 



Using Eq. ( |9b|) and Eq. ( 1T61) . we found 
f(z,t) = f(z,t) + d z - 



{9i + 92) f + ^9i9i ■ d z f 



Collecting Eq. flTHJ and Eq. flU}, we found 

v rt\ n f\ -(j9i+92) 1 9i9i 1 , 

By introducing the nonlinear turbulence scattering term T(f), the kinetic equation can 



(18a) 
(18b) 

(19) 
(20) 
(21) 



be written as 



V (f)=T(f) = -d M 



(22) 
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The phase space convection vector a-y and the symmetric diffusion tensor dq- are given 



by 



-9i 



L T.l 



l T,2 



+ 



-92 



At At 



9x9i 
At 



(23) 



(24) 



The finiteness of dj- and a-j-,2 m T{f) is due to the stochasticity developed by the nonlinearity 
on the time scale At; similarly, the finiteness of dc in C(f) is due to the randomness of the 
collisional scattering events on the time scale longer than the typical correlation time in 
collisional process At this point, the scale of At should be further clarified in connection 
with the nonlinear stochasticity 6|. Let 5t be the turbulence correlation time. Note the fact 
that the stochasticity displays on a time scale longer than St, while in a time scale shorter 
than St, the stochasticity does not display itself. If one looks on a time scale At > St, the 
randomness due to the stochasticity makes dj- and aq- 2 finite; if one looks on the time scale 
At <C St, df and aj-,2 are indeed unimportant. Clearly, dj- and a-j-,2 describe the turbulent 
dissipation on the time scale longer than St. In fact, if one determines to resolve the time 
scale much shorter than St, dj- and a-j-,2 can be ignored, and the 0,7-1 ■ d z f term retained 
in Eq. ( 1221) describes all the nonlinearity; this is essentially the underlying physics behind 



Frieman- Chen's equation 



10[ in gyro-averaged form; note that on this time scale, clt,\ is the 



first order velocity in phase space. The robustness of the turbulence scattering term shall 
be demonstrated in the following applications. 

Note that the temporal accumulation of g\ (the infinitesimal I-transform) is G\ (the 
long-time I-transform). It is not hard to show that the present theory recovers the linear 
theory and quasilinear theory. To do this we separate the distribution function into the 
fluctuation part / and the ensemble averaged part (/), that is f — (f) + f, (f) = 0. Note 
that (a.7-,1) = 0) an d (/) evolves on the transport time scale. 

The linear version of the theory is given by 



V / 



-9i 
At 



(f) 



Its solution is clearly 



h = G 1 -d z (f). 
This recovers the previous linear theory [9], [22I . 



(25) 



(26) 



s 



The quasilinear version is given by 



dt(f) 



(27) 



Substituting Eq. (126|) . one finds the quasilinear transport equation 

1 



d t (f) + d z 



■MGid) ■ d z (f) 



0. 



(28) 



This recovers the previous quasilinear theory 7|, |9J . 
The nonlinear transport equation can be given by 



W) + d z 



(a T ,2) - \{d T ) ■ d z 



+ d z 



(a T J) 



(29) 



where the last term reduces to the quasilinear transport in the quasilinear regime according 
to Eq. ([27D and Eq. (|28|). 

The nonlinear diffusion tensor (df) and convection vector (<xr,2) are attributed to the 
stochastic turbulence scattering effects in the nonlinear regime; they are important only 
for the stochastic motion. Comparing Eq. (|26|) and Eq. (1281) . one concludes that the 
quasilinear transport is due to the long-time phase coherence of the particle excursion G\ 
and the distribution fluctuation /. In the nonlinear regime, this phase coherence is broken by 
the developed stochasticity. Therefore, one expects that the stochastic transport dominates 
in the nonlinear regime, and the quasilinear transport dominates in the quasilinear regime. 

The symmetry of the phase space quasilinear diffusion tensor has been used to demon- 
strate the Onsager symmetry relation of quasilinear transport js]] , as is similar to the Onsager 
relation of collisional transport due to the symmetry of collisional velocity space diffusion 
tensor j^, Q]. Noting the symmetry of the turbulent diffusion tensor in phase space, one 
expects that the Onsager relation still holds in the nonlinear regime, an important point for 
experimental data analysis. 

The nonlinear evolution of the fluctuating part of distribution is given by 



[Vo + {a r> 2)-d z ]f = -d z 



a T Af)-\{d r )-d z )f 



+ 



(30) 



Note that a turbulent thermal conduction term [similar to (dj-)] was heuristically added to 
the drift kinetic equation [similar to Eq. (I30|) ] by Kadomtsev and Pogutse to solve for 
the saturation level of /, by invoking Landau's dissipation cascade argument. 
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For simplicity, we introduce the eikonal ansatz, / = fj t e~ t ^ ult ~ k ' x \ When the turbulence 
saturated, Eq. (|30|) should balance. The real terms are from the right-hand side; the first 
term ~ jL,kf is the linear driving term, wit h jr k the linear growth rate; and the seconc. 
term ~ — (d^)k 2 f is the nonlinear damping 23] or the nonlinear resonance broadening jl] 
due to the turbulent diffusion. Balancing the real terms results to the well-known formula 



0E 



23] of turbulent diffusion coefficient 

Kdr>l 



lL,k 

k 2 ' 



(31) 



which can also be obtained by invoking Prandtl's "mixing- length" argument Q, Q]- The 
discussions presented here and following Eq. f)24p make it clear that the new turbulence 
scattering term collects the mysterious dissipation term in Kadomtsev-Pogutse's equation 
to Frieman-Chen's equation. 

Balancing the imaginary terms of Eq. (130]) results to the nonlinear resonance condition 



k ■ v — k ■ 



(a T ,2) ~ ■ ( d r) 



0. 



(32) 



This indicates a nonlinear phase shift effect, which can be related to important topics widely 
discussed in nonlinear optics 24( and plasma physics 25[: the nonlinear frequency shift (e. 
g. chirping) of a wave resonantly excited with a fixed k, and the nonlinear k spectrum shift 
(e. g. self focusing) in propagation of a wave excited with a fixed u. This nonlinear phase 
shift may be related to the ponderomotive effects {J. 

Although the above important points, from the linear and quasilinear to the nonlinear 
physics, can be successfully identified from the turbulence scattering term, it is still a com- 
plicated nonlinear term; however, it is similar to the collisional scattering term. With the 
two scattering terms and a possible driving source term S included, the kinetic equation can 
be written as 

V (/)=T(/)+C(/) + S. (33) 

It is useful to point out that the new method involves only integrating along the unper- 
turbed orbit, instead of integrating along the full orbit in the conventional method. This may 
bring about advantages by making use of our knowledge of the well-understood unperturbed 
orbit, especially for the toroidal axisymmetric tokamak fusion system 22j, I26M28] . It is not 
hard to extend this theory to the gyrokinetic version that decouples the fast gyro-motion 



for the low-frequency turbulence [9 



12|. 



10 



In conclusion, the nonlinear turbulence scattering term has been found by using the I- 
transform method to decouple the turbulent motion from the unperturbed motion. A unified 
kinetic theory has been formulated to include both the turbulence scattering effect and the 
collisional scattering effect on plasma transport. The evolution of the plasma distribution 
function is governed by the ballistic motion of charged particles along the orbit determined 
by the averaged electromagnetic fields, with the effects of fluctuating part of fields described 
by the turbulence scattering term, from which various important physics scattered in the 
literature can be identified in a unified way. Particularly, the proposed turbulence scattering 
term reveals the collections between two widely used equations: Kadomtsev-Pogutse's equa- 
tion with a mysterious dissipation term and Frieman-Chen's equation without the explicit 
dissipation term. 

This work was supported by the National Natural Science Foundation of China under 
Grant No. 11175178. 
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